Search for: All records

Abstract Once per ≈10⁴–10⁵yr, an unlucky star may experience a close encounter with a supermassive black hole (SMBH), partially or fully tearing apart the star in an exceedingly brief, bright interaction called a tidal disruption event (TDE). Remnants of partial TDEs are expected to be plentiful in our Galactic center, where at least six unexplained, diffuse, star-like “G objects” have already been detected, which may have formed via interactions between stars and the SMBH. Using numerical simulations, this work aims to identify the characteristics of TDE remnants. We take 3D hydrodynamic FLASH models of partially disrupted stars and map them into the 1D stellar evolution code MESA to examine the properties of these remnants from tens to billions of years after the TDE. The remnants initially exhibit a brief, highly luminous phase, followed by an extended cooling period as they return to stable hydrogen burning. During the initial stage (≲10⁵yr) their luminosities increase by orders of magnitude, making them intriguing candidates to explain a fraction of the mysterious G objects. Notably, mild TDEs are the most common, and result in the brightest remnants during this initial phase. However, most remnants exist in a long-lived stage where they are only modestly offset in temperature and luminosity compared to main-sequence stars of equivalent mass. Nonetheless, our results indicate remnants will sustain abnormal, metal-enriched envelopes that may be discernible through spectroscopic analysis. Identifying TDE survivors within the Milky Way could further illuminate some of the most gravitationally intense encounters in the Universe. 

Abstract The detection of GW170817/AT2017gfo inaugurated an era of multimessenger astrophysics, in which gravitational-wave and multiwavelength photon observations complement one another to provide unique insight into astrophysical systems. A broad theoretical consensus exists, in which the photon phenomenology of neutron star mergers largely rests upon the evolution of the small amount of matter left on bound orbits around the black hole or massive neutron star remaining after the merger. Because this accretion disk is far from inflow equilibrium, its subsequent evolution depends very strongly on its initial state, yet very little is known about how this state is determined. Using both snapshot and tracer particle data from a numerical relativity/MHD simulation of an equal-mass neutron star merger that collapses to a black hole, we show how gravitational forces arising in a nonaxisymmetric, dynamical spacetime supplement hydrodynamical effects in shaping the initial structure of the bound debris disk. The work done by hydrodynamical forces is ∼10 times greater than that due to time-dependent gravity. Although gravitational torques prior to remnant relaxation are an order of magnitude larger than hydrodynamical torques, their intrinsic sign symmetry leads to strong cancellation; as a result, hydrodynamical and gravitational torques have a comparable effect. We also show that the debris disk’s initial specific angular momentum distribution is sharply peaked at roughly the specific angular momentum of the merged neutron star’s outer layers, a fewr_gc, and identify the regulating mechanism. 

Abstract Triple stellar systems allow us to study stellar processes that cannot be attained in binary stars. The evolutionary phases in which the stellar members undergo mass exchanges can alter the hierarchical layout of these systems. Yet, the lack of a self-consistent treatment of common-envelope (CE) in triple-star systems hinders the comprehensive understanding of their long-term fate. This paper examines the conditions predicted around binaries embedded within CEs using local 3D hydrodynamical simulations. We explore varying the initial binary separation, the flow Mach number, and the background stellar density gradients as informed by a wide array of CE conditions, including those invoked to explain the formation of the triple system hosting PSR J0337+1715. We find that the stellar density gradient governs the gaseous drag force, which determines the final configuration of the embedded binary. We observe a comparable net drag force on the center of mass but an overall reduction in the accretion rate of the binary compared to the single-object case. We find that, for most CE conditions, and in contrast to the uniform background density case, the binary orbital separation increases with time, softening the binary and preventing it from subsequently merging. We conclude that binaries spiraling within CEs become more vulnerable to disruption by tidal interactions. This can have profound implications on the final outcomes of triple-star systems. 

Abstract Disk-fed accretion onto neutron stars can power a wide range of astrophysical sources ranging from X-ray binaries, to accretion-powered millisecond pulsars, ultraluminous X-ray sources, and gamma-ray bursts. A crucial parameter controlling the gas–magnetosphere interaction is the strength of the stellar dipole. In addition, coherent X-ray pulsations in many neutron star systems indicate that the star's dipole moment is oblique relative to its rotation axis. Therefore, it is critical to systematically explore the 2D parameter space of the star's magnetic field strength and obliquity, which is what this work does, for the first time, in the framework of 3D general-relativistic magnetohydrodynamics. If the accretion disk carries its own vertical magnetic field, this introduces an additional factor: the relative polarity of the disk and stellar magnetic fields. We find that depending on the strength of the stellar dipole and the star–disk relative polarity, the neutron star's jet power can either increase or decrease with increasing obliquity. For weak dipole strength (equivalently, high accretion rate), the parallel polarity results in a positive correlation between jet power and obliquity, whereas the antiparallel orientation displays the opposite trend. For stronger dipoles, the relative-polarity effect disappears, and jet power always decreases with increasing obliquity. The influence of the relative polarity gradually disappears as obliquity increases. Highly oblique pulsars tend to have an increased magnetospheric radius, a lower mass accretion rate, and enter the propeller regime at lower magnetic moments than aligned stars. 

Abstract Many studies have found that neutron star mergers leave a fraction of the stars’ mass in bound orbits surrounding the resulting massive neutron star or black hole. This mass is a site ofr-process nucleosynthesis and can generate a wind that contributes to a kilonova. However, comparatively little is known about the dynamics determining its mass or initial structure. Here we begin to investigate these questions, starting with the origin of the disk mass. Using tracer particle as well as discretized fluid data from numerical simulations, we identify where in the neutron stars the debris came from, the paths it takes in order to escape from the neutron stars’ interiors, and the times and locations at which its orbital properties diverge from those of neighboring fluid elements that end up remaining in the merged neutron star. 

Abstract A black hole (BH) traveling through a uniform, gaseous medium is described by Bondi–Hoyle–Lyttleton (BHL) accretion. If the medium is magnetized, then the black hole can produce relativistic outflows. We performed the first 3D, general-relativistic magnetohydrodynamic simulations of BHL accretion onto rapidly rotating black holes using theH-AMRcode, where we mainly varied the strength of a background magnetic field that threads the medium. We found that the ensuing accretion continuously drags the magnetic flux to the BH, which accumulates near the event horizon until it becomes dynamically important. Depending on the strength of the background magnetic field, the BHs can sometimes launch relativistic jets with high enough power to drill out of the inner accretion flow, become bent by the headwind, and escape to large distances. For stronger background magnetic fields, the jets are continuously powered, while at weaker field strengths they are intermittent, turning on and off depending on the fluctuating gas and magnetic flux distributions near the event horizon. We find that our jets reach extremely high efficiencies of ∼100%–300%, even in the absence of an accretion disk. We also calculated the drag forces exerted by the gas onto to the BH and found that the presence of magnetic fields causes the drag forces to be much less efficient than in unmagnetized BHL accretion. They can even sometimes become negative, accelerating the BH rather than slowing it down. Our results extend classical BHL accretion to rotating BHs moving through magnetized media, and demonstrate that accretion and drag are significantly altered in this environment. 

Abstract Short γ -ray burst (sGRB) jets form in the aftermath of a neutron star merger, drill through disk winds and dynamical ejecta, and extend over four to five orders of magnitude in distance before breaking out of the ejecta. We present the first 3D general-relativistic magnetohydrodynamic sGRB simulations to span this enormous scale separation. They feature three possible outcomes: jet+cocoon, cocoon, and neither. Typical sGRB jets break out of the dynamical ejecta if (i) the bound ejecta’s isotropic equivalent mass along the pole at the time of the BH formation is ≲10 −4 M ⊙ , setting a limit on the delay time between the merger and BH formation, otherwise, the jets perish inside the ejecta and leave the jet-inflated cocoon to power a low-luminosity sGRB; (ii) the postmerger remnant disk contains a strong large-scale vertical magnetic field, ≳10 15 G; and (iii) if the jets are weak (≲10 50 erg), the ejecta’s isotropic equivalent mass along the pole must be small (≲10 −2 M ⊙ ). Generally, the jet structure is shaped by the early interaction with disk winds rather than the dynamical ejecta. As long as our jets break out of the ejecta, they retain a significant magnetization (≲1), suggesting that magnetic reconnection is a fundamental property of sGRB emission. The angular structure of the outflow isotropic equivalent energy after breakout consistently features a flat core followed by a steep power-law distribution (slope ≳3), similar to hydrodynamic jets. In the cocoon-only outcome, the dynamical ejecta broadens the outflow angular distribution and flattens it (slope ∼1.5). 

The engulfment of substellar bodies (SBs) such as brown dwarfs and planets has been invoked as a possible explanation for the presence of SBs orbiting subdwarfs and white dwarfs, rapidly rotating giants, and lithium-rich giants. We perform three-dimensional hydrodynamical simulations of the flow in the vicinity of an SB engulfed in a stellar envelope. We model the SB as a rigid body with a reflective boundary because it cannot accrete. This reflective boundary changes the flow morphology to resemble that of engulfed compact objects with outflows. We measure the drag coefficients for the ram pressure and gravitational drag forces acting on the SB, and use them to integrate its trajectory during engulfment. We find that SB engulfment can increase the stellar luminosity of a 1M⊙ star by up to a few orders of magnitude for timescales of up to a few thousand years when the star is ≈10R⊙ and up to a few decades at the tip of the red giant branch. We find that no SBs can eject the envelope of a 1M⊙ star before it evolves to ≈10R⊙ . In contrast, SBs as small as ≈10MJup can eject the envelope at the tip of the red giant branch, shrinking their orbits by several orders of magnitude in the process. The numerical framework we introduce here can be used to study the dynamics of planetary engulfment in a simplified setting that captures the physics of the flow at the scale of the SB. 

The engulfment of substellar bodies (SBs) such as brown dwarfs and planets has been invoked as a possible explanation for the presence of SBs orbiting subdwarfs and white dwarfs, rapidly rotating giants, and lithium-rich giants. We perform three-dimensional hydrodynamical simulations of the flow in the vicinity of an SB engulfed in a stellar envelope. We model the SB as a rigid body with a reflective boundary because it cannot accrete. This reflective boundary changes the flow morphology to resemble that of engulfed compact objects with outflows. We measure the drag coefficients for the ram pressure and gravitational drag forces acting on the SB, and use them to integrate its trajectory during engulfment. We find that SB engulfment can increase the stellar luminosity of a 1M⊙ star by up to a few orders of magnitude for timescales of up to a few thousand years when the star is ≈10R⊙ and up to a few decades at the tip of the red giant branch. We find that no SBs can eject the envelope of a 1M⊙ star before it evolves to ≈10R⊙ . In contrast, SBs as small as ≈10MJup can eject the envelope at the tip of the red giant branch, shrinking their orbits by several orders of magnitude in the process. The numerical framework we introduce here can be used to study the dynamics of planetary engulfment in a simplified setting that captures the physics of the flow at the scale of the SB. 

ABSTRACT The discovery of GRB 170817A, the first unambiguous off-axis short gamma-ray burst (sGRB) arising from a neutron star merger, has challenged our understanding of the angular structure of relativistic jets. Studies of the jet propagation usually assume that the jet is ejected from the central engine with a top-hat structure and its final structure, which determines the observed light curve and spectra, is primarily regulated by the interaction with the nearby environment. However, jets are expected to be produced with a structure that is more complex than a simple top-hat, as shown by global accretion simulations. We present numerical simulations of sGRBs launched with a wide range of initial structures, durations, and luminosities. We follow the jet interaction with the merger remnant wind and compute its final structure at distances ≳1011 cm from the central engine. We show that the final jet structure, as well as the resulting afterglow emission, depends strongly on the initial structure of the jet, its luminosity, and duration. While the initial structure of the jet is preserved for long-lasting sGRBs, it is strongly modified for jets barely making their way through the wind. This illustrates the importance of combining the results of global simulations with propagation studies in order to better predict the expected afterglow signatures from neutron star mergers. Structured jets provide a reasonable description of the GRB 170817A afterglow emission with an off-axis angle θobs ≈ 22.5°.